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We develop an approach for calculating stationary states of rotating Bose-Einstein condensates 
in harmonic traps which is applicable for arbitrary ratios of the rotation frequency to the transverse 
frequency of the trap uj±. Assuming the number of vortices to be large, we write the condensate 
wave function as the product of a function that describes the structure of individual vortices times 
an envelope function, varying slowly on the scale of the vortex spacing. By minimizing the energy, 
we derive Gross-Pitaevskii equations that determine the properties of individual vortices and the 
global structure of the cloud. For low rotation rates, the structure of a vortex is that of an isolated 
vortex in a uniform medium, while for rotation rates approaching the frequency of the trap (the 
mean field quantum Hall regime), the structure is that of the lowest p-wave state of a particle in a 
harmonic trap with frequency u>±. The global structure of the cloud is determined by minimizing 
the energy with respect to variations of the envelope function; for conditions appropriate to most 
experimental investigations to date, we predict that the transverse density profile of the cloud will 
be of the Thomas- Fermi form, rather than the Gaussian structure predicted on the assumption that 
the wave function consists only of components in the lowest Landau level. 



I. INTRODUCTION 

Bose-Einstein-condensed atomic gases are very well suited to investigating quantized vortex lines. Single vortex 
lines were first made in atomic condensates by Matthews et al. [1], who induced rotation by phase imprinting in a 
spinor condensate. Subsequently, arrays containing many vortices were created in scalar condensates by inducing 
rotation mechanically, either by stirring the condensate [2,3], or by evaporating particles [4]. For a theoretical review, 
see Ref. [5]. In a seminal work, Ho [6] predicted that clouds of atoms confined in harmonic traps, when rotated at 
frequencies close to the transverse frequency lj±_ of the trap, should condense into the lowest Landau level (LLL) in the 
Coriolis force, similar to charged particles in the quantum Hall regime. This insight has led to extensive experimental 
studies in which rotation rates in excess of 0.99 lo± have been achieved, and the structure of the condensate within a 
single cell of the vortex lattice has been examined [4,7]. 

To date, most theoretical work on vortices in harmonically trapped condensates rotating at frequencies close to u>± 
has been based on the use of wave functions in which particles occupy only the lowest Landau level. In contrast, 
for slowly rotating condensates, the usual approach to calculating vortex structure is to solve the Gross-Pitaevskii 
equation. In this paper we address the question of how this approach goes over to the mean field quantum Hall 
description when the rotation rate is increased [8]. We develop a unified method for calculating both the structure 
of individual vortices and the global structure of the cloud for arbitrary rotation rates. Writing the condensate wave 
function as a product of a slowly varying envelope function that determines the density averaged over a single cell 
of the vortex lattice, and a function that determines the variations of the wave function on length scales of order 
the vortex separation and core size, we derive, in Sec. II, the energy of the system. Then in Sec. Ill we derive the 
equation for the structure of the wave function within a single cell of the vortex lattice by variation of the energy 
functional. In Sec. IV we derive equations for the global structure of the cloud. We find that a wave function made 
up only of components belonging to the lowest Landau level is a good approximation only under a very restricted 
range of conditions; if in the non-rotating system the density profile in the plane transverse to the rotation axis is 
of the Thomas-Fermi form, an inverted parabola, then at high rotation the shape remains Thomas-Fermi, and thus 
includes components from higher Landau levels, rather than the Gaussian shape predicted by the LLL calculation. 
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II. BASIC FORMALISM 



We consider a system of weakly interacting bosons trapped in a harmonic potential, V(r) — \m(iJ 2 r\ + ui z z 2 ), 
where rj_ = (x,y), rotating at angular velocity f2 about the z axis. The angular momentum of the system is carried 
in quantized vortices, of number N v » 1 at large rotation rates. We assume the vortices to be rectilinear and to 
form a triangular lattice. When N v » 1 the rotation rate of the lattice is related to the (two dimensional) density of 
vortices, n v , by 

mil = nhn v , (1) 

as follows from identifying the mean velocity via the condition for quantization of vorticity, § c d£ ■ v — hN v (C)/m, 
where N V (C) is the number of vortices surrounded by the contour C. We will not address here the form of the 
corrections to Eq. (1) for finite N v [9]. 

In order to separate out the short distance vortex structure from the large scale structure, we follow the approach 
of Fischer and Baym [10] and write the order parameter as 

V,(r ) = e^f(f )y/W), (2) 

the product of a rapidly varying real factor, f(r), which vanishes at each vortex core, times a slowly varying real 
envelope function, y/n(f), and a phase factor. We normalize f 2 so that it averages to unity over each unit cell of the 
lattice; thus n(r ) is the smoothed density profile of the system, which varies slowly over the unit cells of the vortex 
lattice. The factor e 1 */ describes the local swirling of the fluid - with the phase $ wrapping by 2tt around each vortex 
- together with the overall rotation of the vortex lattice at f2. We generally set K = 1. 
The total energy of the system in the laboratory frame is 

E = |d 3 r|^-|V^| 2 + T/(rMr)/ 2 (r) + |n 2 (r)/ 4 (r)|, (3) 

where we assume a two body interaction described by an s-wave scattering length, a s , with g = ATia s T% 2 /m. With 
Eq. (2), the kinetic energy in the laboratory frame becomes 

/ d3r i |v,A|2 = K = J d3r i {( v ^) 2 + ( v$ ) v 2 + n( yff + \^f 2 • v ")} • ( 4 ) 

We integrate the final term by parts to give — | J f 2 V 2 n\ since n varies slowly across a unit cell of the vortex lattice, 
we may replace the / 2 here by its average (=1) in the cell, so that the integral gives only a vanishing surface term. 
Thus 

K = J d\±- {(VV^) 2 + (V$) 2 n/ 2 + n(V/) 2 } . (5) 

In the unit cell centered on vortex j at position Rj in the plane transverse to the rotation axis, the velocity V$/m 
is the sum of the solid body rotation of the position of the vortex, fix Rj, plus the local velocity around the vortex, 
which we write as V(f>j/m: 

V$(r) ~ infix Rj + W(j)j. (6) 

The (V$) 2 term thus becomes 

/ ^ (V#)» - Y^jfmf \SM- + I m!! ^} ; (7) 

the integration is over unit cell j, and the sum is over all cells. The cross term vanishes since in the limit N v 3> 1 the 
average velocity in the cell measured with respect to the center of the cell vanishes. In the Wigner-Seitz approximation, 
which we employ below, <j>j becomes the azimuthal angle measured with respect to the point Rj . Writing within cell 
j, Rj = r± — p, the final term in Eq. (7) becomes \lil 2 — J2j Ij jinVi 2 p 2 nf 2 , where I = J d?rmnf 2 r\ is the total 
moment of inertia of the system. Similarly the transverse trapping potential term becomes \lu 2 ■ 
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To determine the equilibrium structure, we work in the frame rotating at angular velocity f2. [This procedure is 
equivalent to determining the equilibrium structure at fixed angular momentum, L, by minimizing the total energy 
taking the constraint of fixed L into account by a Lagrange multiplier, ft.] The total angular momentum along the z 
axis is given by 



L = j d 3 rn(r)f 2 (r) (fx V*(r)) z . 
Using Eq. (6), and again writing in cell j, f± = Rj + p, we have, 

L = m + J2 J d 3 rnf{{px V<t>j) z - rnQp 2 } 



(8) 



(9) 



since the average position in the unit cell j is Rj, and the average velocity in the local frame of the vortex j vanishes. 
The first term is the angular momentum of the center of mass of the cell, and the second the intrinsic angular 
momentum within the cell. Assembling the pieces, Eqs. (3), (5), (7), and (9), we have, 



E' = E-QL = J 



d A r 



V 2m 



rn 



[ d 3 rn 



(10) 



Expressing I in terms of the moment of inertia, I — J d 3 rmnr 2 _, of the smoothed density distribution we write 

E ' = S d * r {^r + ? n(r)w * v ) + ° 2 ) J+ E^- (") 



where 



Ei = \f rn {^Sr + £ (v ^ )2 + ^p 2 ^ 2 + \ mn2 p 2 ^ x v< ^ + 9 2 nf4 

is the internal energy within cell j. 



(12) 



III. EQUILIBRIUM STRUCTURE OF VORTICES 

We turn now to determining the structure of the vortices within the unit cells. To do this we introduce the Wigner- 
Seitz approximation to evaluate the vortex sum, replacing the hexagonal unit cell by a circle of radius £ = l/(mf2) 1 / 2 . 
Then / is cylindrically symmetric within each cell. In the following, we assume that the vortex spacing is small 
compared with the characteristic length scale in the axial direction. The term in Eq. (12) containing df/dz can 
then be neglected, and / depends only on the transverse coordinate and the average local density. In cell j, <j>j 
becomes the azimuthal angle with respect to the center of the cell. Again we write within cell j, f± = Rj + p, so 
that (V</>j) 2 = l/p 2 . Furthermore, (p x V<t>j) z becomes just K, so that the angular momentum in the Wigner-Seitz 
approximation is, 

L = m + J2 J^d 3 rnf 2 (l - mflp 2 ) = Ifl + N - , (13) 

where (p 2 ) ~ £ 2 /2 is the average of p 2 within a given cell. We have neglected gradients of the smoothed density 
here. For an incompressible fluid, (p 2 }/£ 2 = 1/2, and therefore the additional angular momentum per particle in 
the Wigner-Seitz approximation is h/2, which is close to Tkachenko's result [11,12] for a triangular lattice in an 
incompressible fluid, (7r/4\/3)7i ~ 0.453?i. 

In the Wigner-Seitz approximation, Ej becomes, 

+ ^lp 2 (f 2 1) + \mtfp 2 -n+ 9 -nA . (14) 



En 



f-L 






| 2m 




p 2 
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The form of / within each cell is determined by minimizing Ej with respect to /, subject to J\ cPpf 2 — 1, with the 
boundary conditions that /(0) = and df /dp — at p — £. Since there are no terms coupling / at different values 
of z, the equilibrium / depends on z only through the dependence of the average density on z. Thus within a given 
cell, at given height z, 

M~jh + ?) + + snf ' »-<"«"'»'- <15) 

This equation describes the vortex structure for all values of parameters, provided that N v is large. Equations (10) and 
(15) generalize the result of Ref. [10] through inclusion of the (\7y/n) 2 /2m and mus\p 2 f 2 /2 terms. The (V-^) 2 term 
allows us to go beyond Thomas- Fermi, when this energy dominates the interaction term. In the limit fl 3> u)j_/2gn, 
appropriate to the regime described in Ref. [10], the u>\ term in (15) can be neglected. 
It is useful to define the averages over the unit cell, 




^=2^<P 2 (/ 2 -l)> 



(17) 



and 

b={f% 

these quantities are dependent on the density within the cell. Then quite generally, 



E, 



(18) 



(19) 



For slow rotation, the core structure is basically that of a single vortex [13], and is reasonably well approximated 
by [14] 



mi+P 2 ) 1 ' 2 ' 



(20) 



where £o = h/ ^/2mgn is the Gross-Pitaevskii healing length. The corresponding density within this approximation 
to / is shown as curve a in Fig. 1 for the particular value, £o = 0.1£. 




0.2 0.4 0.6 

Radius in units of cell radius 



FIG. 1. Density within a vortex core in units of the average density in the cell, as a function of the transverse radius in units 
of the core radius I: a) the single vortex form, (20); b) the linear core approximation, (21); and c) the quantum Hall structure, 
(29), d) the free particle Bessel function Ji (dashed). Curves a and b are calculated for £o = 0.1£. 
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Reference [10] used a simple linear approximation for / for all rotation speeds, in which / rises linearly to the 
effective core radius £ and then becomes constant to the edge of the cell, 



fip) 



1 



(1 -£2/2^ 2 ) 1/2 



i, t<p<t 



(21) 



The corresponding density is shown as curve b in Fig. 1, for the value £ = \/6£o with £o = 0.1£ In general one can 
solve Eq. (15) numerically for s, although we shall not do this here. 

With the linear approximation (21), the individual vortex energy, Ej, becomes 



E, 



<l'm[na(0 + ^a h (0 + ^b(0-^ 



where 



a(C) 



l_jdnC 

1 - C/2 ' 



3C ~ 2C 2 
12(2-C)' 



(22) 



(23) 



and C = £ 2 /^ 2 is the fractional area occupied by the vortex core. The fluctuations in the density within a cell 
renormalize the (long wavelength) coupling constant [10] by a factor b = (n 2 ) / (n) 2 > 1, given, for the ansatz (21), by 



6(0 



1 - 2C/3 
(1 - C/2) 2 



(24) 



The relative area occupied by the core at position (rj_, z) is found by minimizing the integrand of (22) at the density 
n(r±,z): 



In the Thomas-Fermi regime, the sound velocity, s, in the center of the trap is given by 

2/5 



ms 2 = gbn(0) = 



15Nba, 



d± uj± 



J- r> 



(25) 



(26) 



where d± = l/frnwi) 1 / 2 is the oscillator length for transverse motion. We show, in Fig. 2, the corresponding 
prediction for £ at the center of the trap as a function of rotational velocity for 87 Rb, taking the representative values, 
N = 2.5 x 10 6 , w_l/2tt = 8.3 Hz, and w z /2tt = 5.2 Hz. 




FIG. 2. Variation of the core size with rotational velocity, in the linear approximation to the core structure. 
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As we see in Fig. 2, the core structure changes rapidly as f2 approaches the transverse trap frequency u±. In order 
to study rotational velocities comparable to the transverse trapping frequency it is useful to spread out the horizontal 
scale by measuring rotational rates in terms of the rotational rapidity, y, defined by [15]: 



tanhy, 



(27) 



1, wi + 
y= 2 ln ^fi- 



(28) 



The rapidity variable essentially counts the number of 9's in the fraction f2/wj_ as the fraction approaches unity (just 
as metal dealers describe the purity of metals). For example, the currently achieved [7] Q/u)j_ = 0.995 corresponds to 
a rapidity of 3.00, while Q/uj± = 0.999 corresponds to y — 3.45, and fl/u}± = 0.9999 to y = 4.61. In Fig. 3 we show 
the variation of the core size in Fig. 2, now as a function of rapidity. 
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4 6 

Rotational Rapidity 

FIG. 3. Variation of the core size shown in Fig. 2, as a function of rotational rapidity, y. 



In the mean field quantum Hall regime, where Q approaches ui±, the cloud expands to the point where the gn 
term becomes a small perturbation on the structure within a cell; to lowest order / assumes the particularly simple 
oscillator p-state structure, 

/ = C ^ e -p 2 /2^ (29) 

plus small terms, where C = (1 — 2/e) -1 / 2 . The p-wave solution for a particle in a trap with frequency f2 which 
satisfies the usual boundary condition / — * for p — ► oo kindly has zero slope precisely at p = I. With this form of /, 



Ei =jd 3 m (SI + - ^V(/ 2 - 1) + f nf 



(30) 



so that 



this energy becomes the exact quantum Hall result in the limit, fl — > ujj_. The average of p 2 /£ 2 is given by 

in the linear core approximation one finds instead, 0.614. The renormalization of the coupling constant by fluctuations 
in the density within a cell, the factor b = (n 2 )/(n) 2 , is given in the quantum Hall limit by [16,17], 
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6-1^-1-^1 - 1 158- (33) 

by comparison, the linear core approximation yields 1.192. Note that the total moment of inertia, /, equals I + N(b' — 
l/2)/fi. 

IV. GLOBAL STRUCTURE OF THE SYSTEM 

We turn now to determining the global structure of the cloud, by minimizing 

E'-Jd'rf^^ + ^A+^l-^I + ^E,, (34) 

3 

at fixed particle number. From Eq. (34) we derive the effective Gross-Pitaevskii equation for the smoothed density: 

{V 2 m 1 
— + — (co 2 z z 2 + {uj\ - n 2 )rl) + M C eii(n) J. = liy/n, (35) 

where /i cc n = SEj/Sn, defined by Eq. (15), contains the explicit interaction energy term gnb. In the quantum Hall 
limit, Eq. (30) implies that 

,lo 2 -n 2 

Mceii = 20 + bgn + b' x . (36) 

The structure in the axial direction will, for O sufficiently close to iv±, always become Gaussian. The criterion 
for the axial structure to be Gaussian is that gn be small compared with the axial oscillator frequency, u) z . Since 
the system density falls indefinitely with increasing ft, this condition will eventually be satisfied. From Eq. (26), the 
criterion is 



fm z b\ 3/2 2 5 l 2 d ± 
lu 2 , \ uj± J l5Na s 



7^ (37) 



ujj_\ 3/2 15Na s 



or in terms of rapidity, 

For very weak interaction the transverse structure is Gaussian in the quantum Hall limit [6]: 

n(r ± , z) = Ka{z) 2 e- r2 /^ 2 N{z), (39) 

where Af(z) is the number of particles per unit length in the axial direction. As we shall see below, such a Gaussian 
describes the system only for Na s <C d z , where d z = l/(mw 2 )'/ 2 is the axial oscillator length. For this Gaussian, 



plus a constant times N. Minimizing with respect to a(z) at fixed Af(z), we find, 



(40) 



a(z) = ,I(U 2^A-(:),J I '{-^^)'\ (41) 



in agreement with Ref. [6]. 

However, if the transverse structure of the non-rotating cloud is Thomas-Fermi, it will remain Thomas-Fermi as 
the cloud is spun up, even to the quantum Hall limit. The criterion for the transverse structure to be Gaussian is 
different than in the axial direction, since the effective transverse oscillator frequency, (ujj_ — O 2 ) 1 / 2 , goes to zero. 
The criterion becomes instead that the interaction energy, gn, be small compared with the transverse kinetic energy: 
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gn ~ gN/ZR^ <C l/2mi?]_, where Z is the axial height and R± the transverse radius. This condition implies that 
Na s /Z be <C 1. Since the total density per unit height, N/Z, increases with increasing fl as the system flattens 
out, the structure in the transverse direction can only be Gaussian if the transverse structure in the non-rotating 
cloud is itself Gaussian. The maximum that N/Z can become is ~ N/d z , where d z is the axial oscillator length. For 
Na s /d z ;$> 1 the structure in the radial direction will be Thomas- Fermi at large Q, even if it is Gaussian at small Q. 

Note that in the quantum Hall limit, even though the interaction energy plays only a perturbative role in determining 
the structure within each cell of the lattice, it is crucial in determining the global structure. In particular, it is 
responsible for inclusion of components from higher Landau levels required to produce a Thomas-Fermi profile. 

The final axial-Gaussian, transverse-Thomas-Fermi structure at high rotation has the form 



n(O=e-» /<, M»(0)-^( W i-fiM]- (^) 



Using J d 3 rn — N, we find 



7T 



3/2 



N = — d z Rln(0), (43) 



where the transverse size, R±, is given by the point where n(r) falls to zero, 

V2 , ,„ / , ,2 \ 1/2 



MdPW =«.(0»-(^) . (44, 



In terms of the total number, N, 



d ttV 8 V d z uo 2 - n 2 



Nb-^-r^ , (45) 



and 



,2 r>2 \ i/ 2 



m = ^rAJ^r^ L ) • (46) 



27r 5 /4 ybd 4 d z a s uj\ 



V. MEASURING THE CORE SIZE 



Several quantitative measures can be used to compare predicted core sizes with experiment, and with theory in the 
quantum Hall regime. The first is simply to compare the slopes of / at the origin. The slope of the order parameter 
in the linear approximation, Eq. (21), is l/£(C(l — C/2)) 1 ^ 2 , which approaches 1.62/^ as Q — > u>±. On the other hand 
the quantum Hall wave function, (29), has slope 1/1(1 — 2/e) 1 / 2 = 1.95/i. The second is to measure the mean square 
radius, r 2 , of the density deficit in the core, defined by, 

2 = Ij d 2 P [f(e? f( P ) 2 ] p 2 

Tc I 3 d 2 p[.f(l) 2 -.f(p) 2 ] • ( ' 

For the quantum Hall wave function, r 2 /£ 2 = ((11/2) — 2e)/(3 — e) ~ 0.225. To compare with the result from the 
linear approximation to /, we note that r 2 /l 2 = £/3, while as f2 approaches lj±, the value of ( is found from the 
minimum of a(£) +af l (C), which is at ( ~ 0.519; thus in this limit, r 2 /£ 2 ~ 0.173. Note that although the initial slope 
of the quantum Hall wave function is larger than that in the linear approximation, the mean square radius of the 
depression is also larger, since the depression in the quantum Hall wave function extends over the entire cell. Both 
measures of the core size in the linear approximation are in reasonable agreement with the exact quantum Hall result, 
given the simplicity of the approximation. 

Experimentally, core properties are investigated after the rotating cloud has expanded. In the JILA experiments, 
the atoms are transferred to a state in which the magnetic forces tend to drive the cloud apart. It is therefore 
necessary to investigate how the vortex-core structure is affected by the transfer to the new state and the subsequent 
expansion of the cloud. Under expansion, the density drops, eventually reaching the point where the interaction 
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energy no longer plays a role in determining the structure within the individual cells. The centrifugal force plays 
no role within a cell. If the potential is adiabatically turned off, allowing the system to expand slowly to the 
point where the interaction within a cell is small compared with the bending energy of the order parameter within 
a cell, or l/2m£ 2 — N v /2mR 2 L 3> gN/R\Z. which is the case when the axial height expands to the point where 
Z ^> 8Tra s N/N v , then the structure within the individual cells, given by Eq. (15), is the Bessel function, C\J\(x r / £'), 
where x — 1.84 is the location of the first maximum of the J\{x), C\ = 2.05 and £' is the cell size in the expanded 
cloud. In fact, the Bessel function solution is always within 0.015 of the quantum Hall solution, and the two solutions 
would be effectively indistinguishable in practice (see Fig. 1). The slope at the origin of the Bessel function is 1.88, 
compared with 1.85 for the quantum Hall solution, while the mean square radius, r 2 /£ 2 , of the depression is 0.231, 
compared with 0.225 for the quantum Hall solution. 

One can distinguish two stages in the evolution of the cloud during release and the subsequent expansion. The first 
is the period when the atoms arc transferred to an untrapped state, and the second is expansion of the cloud in a 
modified trapping potential. The transfer of atoms occurs on a time scale short compared with dynamical times for 
the particles. Therefore the sudden approximation should be good, and changes in the coordinate-space wave function 
during the transfer should be negligible. This implies that both the global structure of the cloud and the structure of 
an individual cell of the vortex lattice are unchanged. After transfer of atoms to the new state, the structure within a 
cell will not correspond to the equilibrium configuration for the particular rotation rate because of the change in the 
trapping potential, which is determined by the instantaneous value of u±. The calculations described in the previous 
paragraph demonstrate that when interaction effects are small, the structure of the condensate wave function within 
a single cell depends only weakly on u±. Therefore, after transfer, the wave function within a cell will be the lowest 
state for the new value of uj±, apart from corrections of order one per cent. Likewise, for rotation rates so small that 
interaction effects dominate, we expect a similar conclusion to hold because the oscillator potential plays little role in 
determining the structure of an individual vortex. 

We now consider the degree to which the vortex cores adjust adiabatically in the expansion. To do this, we compare 
the time scale, r co n, for response of the structure of a cell of a vortex lattice with the expansion time scale, T cxp . When 
the vortex core radius is small compared with the cell radius, the time for adjustments of the core is of order the core 
radius, <~ {mgn)^ 1 / 2 , divided by the sound speed s, or t cc \\ <~ Ti/gn. When the core radius becomes comparable to 
the cell radius, i.e., Ml <; gn, the inverse response time becomes of order the kinetic energy associated with a particle 
confined within a volume of radius £, divided by Ti, or r cc n ~ m£ 2 /h — l/f2. Thus l/r co n is always the larger of gn/% 
and fi. These estimates should apply at all stages in the evolution, provided that n and are the instantaneous 
values of these quantities. We note that if the expansion is purely two-dimensional, a good approximation for the 
recent experiments [7], the density and Q, both scale as l/R^; therefore, the ratio gn/TiQ, remains constant, and the 
core expansion rate always remains gn/h or O. 

There are similarly two regimes for the expansion. At low rotation rates, when the interaction energy per particle, 
gn, is large compared with fkj±, the expansion velocity is determined by the interaction energy of the cloud, and 
is typically of order the sound velocity, so, in the cloud before release (the subscript denotes quantities just prior 
to release). On the other hand, when the typical initial orbital velocity, £IqR±o, exceeds the sound velocity, sq, the 
dominant contribution to the expansion velocity after switching off the trap potential is the orbital motion, and 
therefore the expansion velocity is of order Q R±o- The typical expansion rate, 1/t cxp , is thus always the larger of 
Sq/R± and Q R±o/R±. 

Now let us compare time scales. For low rotation velocities, fig /$ sq/ R±0i w ^ have 



T cc n 1 noR±o 1 R± 

T cxp ms R±o nR± msoR_Lo R±o 



(48) 



where the latter estimate holds for two dimensional expansion. This ratio is initially smaller than unity, implying that 
the cell initially adiabatically adjusts during the expansion, but if the cloud expands to a radius i?^_ /£o, where £o 
is the Gross-Pitaevskii healing length, the condition for adiabaticity will be violated. For intermediate rotation rates, 
Sq/R±o ~^o~5 n o/^, the ratio of times is given by 

Tcell ^ hfl R±p 

T eX p gn R± ' 

This ratio starts at a value less than unity but increases oc R±/R±o as the cloud expands. For the final case of fast 
rotation, flo^gno/h, the ratio is 

T co n ^ O i?J_0 ^ R± 
r OX p flR± R±Q 
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In this case the adiabatic assumption is marginally satisfied initially, and is violated during the subsequent expansion. 
We conclude that one may draw no general conclusions about the development of vortex core structure during 
expansion on the basis of arguments about time scales; more detailed studies are needed. 

It is interesting to note that states made up only of components in the lowest Landau level expand homologously 
when the effects of interaction arc neglected [18]. In this case the structure of a single cell remains invariant, with 
only changes in scale, independent of the transverse length entering the wave function [cf. Eq. (29)]. Even though the 
condition for adiabaticity is violated, the structure of the single cell is precisely what would be predicted assuming 
adiabatic behavior. 



VI. CONCLUSION 



In this paper we have developed a unified framework for describing the structure of rotating Bose-Einstcin con- 
densates containing a large number of vortices. We have derived a Gross-Pitaevskii equation which describes the 
structure of individual vortices and have demonstrated how the mean field quantum Hall state emerges as a simple 
continuation of the structure for small rotation rates. We find that the global density profile of the rotating clouds 
in the transverse direction is generally of the Thomas-Fermi form, rather than the Gaussian that emerges if only the 
lowest Landau level is occupied. 

A number of open problems remain for future work. Throughout, we have assumed that the Gross-Pitaevskii 
approach may be used, and have neglected effects of excited states. One such effect is the zero-point motion of 
collective modes [21,17], which broadens the density profile of individual vortices and makes the lowest density non- 
zero. The density is the center of the vortex can also become non-zero via anomolous modes of excitation of the 
condensate [5] that in the linear approximation have a negative excitation energy. Even at zero temperature, such 
modes will have a non-zero equilibrium population such that the energy of an anomalous modes, including the effects 
of self-interaction, is just equal to zero. Explicit calculations are given in Rcf. [19]. A third effect is the thermal 
population of excited states, which likewise will lead to a non-zero density at the center of the vortex. All of these 
effects must be taken into account in a detailed comparison of experiment with theory. Further problems include 
the quantitative delineation of the effect of expansion on the vortex core structure, and inclusion of effects of the 
lattice beyond the Wigner-Seitz approximation, such as the rigidity to shear motion, which manifests itself, e.g., in 
Tkachenko modes [11,20,21]. 
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